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Abstract 

Bohmian mechanics allows to visualize and understand the quantum-mecha- 
nical behavior of massive particles in terms of trajectories. As shown by 
Bialynicki-Birula, Electromagnetism also admits a hydrodynamical formula- 
tion when the existence of a wave function for photons (properly defined) 
is assumed. This formulation thus provides an alternative interpretation of 
optical phenomena in terms of photon trajectories, whose flow yields a pic- 
torial view of the evolution of the electromagnetic energy density in configu- 
ration space. This trajectory-based theoretical framework is considered here 
to study and analyze the outcome from Young-type diffraction experiments 
within the context of the Arago-Fresnel laws. More specifically, photon tra- 
jectories in the region behind the two slits are obtained in the case where 
the slits are illuminated by a polarized monochromatic plane wave. Expres- 
sions to determine electromagnetic energy flow lines and photon trajectories 
within this scenario are provided, as well as a procedure to compute them in 
the particular case of gratings totally transparent inside the slits and com- 
pletely absorbing outside them. As is shown, the electromagnetic energy 
flow lines obtained allow to monitor at each point of space the behavior of 
the electromagnetic energy flow and, therefore, to evaluate the effects caused 
on it by the presence (right behind each slit) of polarizers with the same or 
different polarization axes. This leads to a trajectory-based picture of the 
Arago-Fresnel laws for the interference of polarized light. 
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1. Introduction 

By means of different experimental techniques based on the use of low- 
intensity beams, nowadays it is possible to observe the gradual emergence 
of interferencepatterns in the case of both massive particles [l|, 0| and light 
(photons) 0, Q, That is, the patterns are built up by detecting single 
events (one by one) which can be followed on a screen as a function of time. 
In the pattern formation sequence, at short time scales, when the number 
of events (or detections) recorded is still relatively low, seemingly random 
distribution of points can be observed on a screen used to monitor the ex- 
periment. However, as time proceeds and the number of detections increases 
sufficiently, the well-known band structure of alternate light and dark inter- 
ference fringes starts to appear. Hence, moving from the darker regions of 
the pattern to the lighter ones means that the density of particles detected 
increases. This kind of experiments constitutes a nice manifestation of the 
statistical nature of Quantum Mechanics, which, in the large count-number 
limit, establishes that both massive particles 0] and photons 0] distribute 
according to a certain probability density. 

Based on this kind of experiments, one could assume that each point 
on the screen corresponds to the final position of the trajectory pursued by 
a detected particle. This immediately leads to considering trajectory-based 
formulations to explain interference processes, where the analysis of the prop- 
erties displayed by (particle) trajectories provides a dynamical explanation 
and interpretation of the process. In other words, these formulations would 
contain the mechanisms that allow to account for the appearance of fringes in 
interference experiments with quantum particles (either massive particles or 
photons) in the large-number limit of single-event counts, this outcome being 
in agreement with the results obtained from standard wave-based formula- 
tions. This is an important and remarkable advantage, for these approaches 
thus constitute physical models where both corpuscular and wave-like prop- 
erties are compatible. 

In the case of interference experiments with massive particles, different 
trajectory-based studies can be found in the literature [g, |9|, M, El, H, 
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14], |15|, ll6[, which are aimed to establish a bridge between the single counts 
observed experimentally and the statistical quantum-mechanical description 
in terms of smooth, continuous probability densities. Many of them rely 



on Bohmian mechanics [17], [18|, [19|, |20|, l2l| . a hydrodynamic formulation of 



Quantum Mechanics in terms of trajectories, which, to some extent, can be 
considered at the same level as Newtonian mechanics with respect to (classi- 
cal) Liouvillian mechanics. Bohmian trajectories evolve following the stream- 
lines associated with the quantum probability current density (quantum flow) 
and, therefore, reproduce exactly the quantum-mechanical probability distri- 
bution in both the near and the far field when a large number of them is 



Alterna- 



considered, in particular, in interference experiments 
tively, the emergence of interference patterns in the far field has also been 
simulated by using trajectories determined from the momentum distribution 
(MD trajectories) associated with the particle wave function [l6| . 

In the case of interference processes with light (photons), similar descrip- 
tions are not so numerous. Some early attempts to explain interference and 
diffraction in terms of electromagnetic energy (EME) flow lines were those of 
Braunbek and Laukien (22| and Prosser 23, 24j . Within this approach, based 
on Maxwell's equations, the photon trajectories are associated with the EME 



flow lines [24J, which are obtained after solving a trajectory equation arising 



from the Poynting energy flow vector [25]. Similarly, but starting from the 



Dirac equation instead of Maxwell's ones, Ghose et al. [26| also determined 
photon trajectories. The formal grounds of these approaches are similar 
to those of the hydrodynamic formulation of Electromagnetism developed 
by Bialynicki-Birula 27], |28|, |29|, |30j. This formalism arises after assuming 



that, as happens with massive particles, photons can also be described by 



a well-defined wave function, this beint 
attention in the literature 27], |28|, |29j, [ 



a problem which has received much 



40, 41, 42 



31|, |32|, |33|, |34|, |35|, |36|, |37j, |38j, [39, 



Specifically, by following the same argumentation that led Dirac 
to the relativistic equation for the electron, Bialynicki-Birula 27, 28, 29, 3ol 
reaches Maxwell's equations. Accordingly, though these equations have a 
purely classical origin, they can be expressed as a Dirac-like equation, with 
the wave function satisfying it being a complex vector encompassing both 



the electric and magnetic field, namely the Riemann-Silberstein [43 
Appendix [A} . 



sec 



More specifically, in order to better understand the relationship between 
EME flow lines and photons, consider the following. As mentioned above, 
the EME flow lines describe the paths along which the EME flows. Assuming 
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that a wave function denned in terms of the EM field can be associated to 
a quantum of energy or photon 27J, |28|, [29 



38. 39 



41, 42 



3o|, i3il, m, l33|, m m \m m, 



the probability that a photon arrives at a certain region 
on a screen will be proportional to the EME density at that region, i.e., to 
the density of flow lines reaching such a region. Hence, one could interpret 
that the trajectory arrivals somehow simulate the arrivals of single photons 
to the screen, as seen experimentally That is, the trajectories are the 

paths along which the EME flows and ensembles of them will constitute the 
eventual/possible paths taken by single photons. However, a photon itself 
has to be represented by an ensemble of trajectories in order to satisfy the 
condition of energy quantization. Throughout this work, we are going to 
use the term "photon trajectory" to denote any of the eventual EME flow 
lines or trajectories that a photon can pursue, but taking into account the 
previous discussion. 



Recently, we have considered [45| Prosser's approach (a simplified view 
of the more general Bialynicki-Birula hydrodynamical formulation) to study 
Young's experiment and the Talbot effect for a monochromatic, linearly po- 
larized incident electromagnetic (EM) field (or wave) reaching the diffracting 
grating. This study was carried out in the region between the grating (be- 
hind it) and the detection screen, thus the photon trajectories were computed 
directly from the diffracted wave function, which was obtained by requiring 
that the components of both the electric and the magnetic vector fields sat- 
isfied Maxwell's equations as well as the boundary conditions imposed by 
the grating. The trajectories obtained showed how the EME redistributes 
in space from the grating to the detection screen, located far away from the 
former in the Fraunhofer region [j^ 10]. By properly sampling the field at 
the slits, one could also observe how, as time proceeds, the accumulation of 
trajectory arrivals, which eventually led to the appearance of the well-known 
fringe-like pattern, just as in the analogous real experiment 0, 0, [Hj or as in 
Quantum Mechanics, where the formation of similar patterns arises after a 
large count of particle trajectories. 

In this work, we go a step forward and present a more general trajectory- 
based analysis of Young-type experiments carried out with polarized light. 
Experiments with high-intensity polarized beams are well-known since the 
beginning of the XlXth century, when Arago and Fresnel enunciated their 
laws for Young's interference experiment with polarized light 46|. Accord- 



ing to these laws, no interference pattern is observed if, for example, the 
two interfering beams are linearly polarized in orthogonal directions 47] or 
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both are elliptically polarized but with opposite handedness [48]. By means 



of the Bialynicki-Birula hydrodynamical formulation of Electromagnetism 



271 . 1281 . 1291 . |30| . here we provide a trajectory picture for this kind of ex- 
periments. Since we are interested in the interference phenomenon and its 
relationship with the polarization of the interfering waves, we start the cal- 
culation of the photon trajectories behind the slits. By assuming that the 
incoming wave is a polarized monochromatic EM plane wave, the components 
of the electric and magnetic fields can be expressed in terms of a function 
that explicitly takes into account the boundary conditions imposed by the 
two slits (here, the grating is considered to be totally transparent inside the 
slits and completely absorbing outside them). Within this scenario, we pro- 
vide expressions to determine the photon trajectories as well as a procedure 
to compute them in the particular case of gratings totally transparent inside 
the slits and completely absorbing outside them. As is shown, the photon 
trajectories obtained allow to monitor at each point of space the behavior of 
the electromagnetic energy flow and, therefore, to evaluate the effects caused 
on it by the presence (right behind each slit) of polarizers with the same or 
different polarization axes. This thus constitutes a very interesting pictorial 
view of the Arago-Fresnel laws for the interference of polarized light, consis- 
tent with the more conventional one based on standard Electromagnetism. 

The organization of this paper is as follows. In order to be self-contained, 
in Section [2] we introduce some fundamental theoretical grounds related to 
the photon-trajectory formalism, starting from Maxwell's equation, as also 
done within the framework of the Bialynicki-Birula hydrodynamical formu- 
lation of Electromagnetism. Particularly, we will focus on those elements 
related to the analysis of experiments with polarized light. In Section [31 
we describe the incident EM field and its polarization, as well as the deriva- 
tion of the corresponding photon trajectories before reaching the interference 
grating. In Section HJ the EM field at and behind the grating is obtained as 
a solution of Maxwell's equations. This field is going to be used to calculate 
later the associated photon trajectories. In Section [5j we present and discuss 
the photon trajectories behind a two-slit grating in the case of an incident 
circularly polarized EM field. In Section an analogous analysis is carried 
out, but when the slits are covered with polarizers with orthogonal polar- 
ization axes. Calculations of photon trajectories for the cases of incident 
linearly and circularly polarized EM fields are shown. Finally, in Section [7] 
the conclusions arisen from this work are summarized. 
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2. Maxwell's equations, EME flow lines and photon trajectories 

Consider a monochromatic EM field propagating in vacuum. Under these 
conditions, the electric and magnetic components of this field can be ex- 
pressed as harmonic waves, as 

E(r, t) = E(r)e^, H(r, t) = H(r)e" ia; ', (1) 

respectively, where the space dependent parts of these fields satisfy the time- 
independent Maxwell's equations, 

V-E(r) = 0, (2) 

V-H(r) = 0, (3) 

V x E(r) = iwjx H(r), (4) 

V x H(r) = -zu;e E(r), (5) 

as well as the boundary conditions associated with the particular problem 
under study. Equivalently, from ©-(jS]), it is readily shown that both E(r) 
and H(r) satisfy the Helmholtz equation, 

V 2 E(r) + k 2 E(r) = 0, (6) 
V 2 H(r) + A; 2 H(r) = 0, (7) 

where k = u/c. 

The EME flow lines are obtained from the real part of the time-averaged 
complex Poynting vector [iil ]. 

S(r) = ~ Re[E(r) x H*(r)], (8) 

as follows. We know that this vector describes the flow of the time-averaged 
EME density through space, 

U(r) = \ [eoE(r) • E*(r) + /i H(r) • BT(r)] . (9) 

That is, the EME density is transported through space as a flow described 
by S(r). Formally, this can be expressed as 

S(r) = E7(r)v, (10) 



6 



where v is a local effective velocity vector field, namely the ray velocity 



49|. This velocity indicates the direction of the EME flow at each space 
point and its magnitude is equal to the EME that crosses in unit time an 
area perpendicular to the flow direction (given by the Poynting vector S(r)) 
divided by the EME per unit volume. Taking into account the Bialynicki- 
Birula hydrodynamical viewpoint, one can assume that the EME travels 
along streamlines defined by ( ITU]) . Therefore, this equation can be recast in 
a more convenient way as 

dr = 1 S(r) 

ds cZ7(r)' 

whose solutions, r(s), are the streamlines or EME flow lines followed by 
the EME density in configuration space or, within a Bohmian-like reinter- 
pretation of the Bialynicki-Birula hydrodynamical formulation, the photon 
trajectories. Here, s is simply a parameter which labels the evolution across 
space of the corresponding trajectory, however, for practical purposes, later 
we will reparametrize the solutions in terms of one of the coordinates — in 
particular, we consider the reparametrization {x,y, z}(s) — » {y, z}(x). 

In the particular case of interference experiments, a functional form for 
pip can be found by means of the following considerations. The screen 
containing the grating is on the XZ plane, at y — 0, and the slits are parallel 
to the z-axis, their width (5) being much larger along this direction than 
along the x-direction (i.e., 5 Z ^> 5 X ). Accordingly, we can assume that the 
EME density is independent of the z-coordinate and, therefore, the electric 
and magnetic fields will not depend either on this coordinate. This allows us 
to consider a simplification in the analytical treatment, for, as mentioned in 



49l | . a problem independent of one Cartesian coordinate is essentially scalar 
and, therefore, can be formulated in terms of one dependent variable. Thus, 
if we introduce here 

§— £ 

into Eqs. (@| and (J5j), we obtain two independent sets of equations, 

dE z iu dE, iuj 

^_ TT ^_ TT 

dy e c 2 x ' dx e c 2 y ' 

dH v dH r „ 

-JT- ~ = -iue E z , 13 

ox oy 
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and 



dH z 8H Z 

— — = -iuie E x , — — = iue E y , 
oy ox 

dE„ dEr iu 



J x 



dx dy eoc 



H z . (14) 



The set (j!3p only involves H x , H y and E ZJ and therefore is commonly referred 
as a case of ^'-polarization, while the set (114)). which only involves E x , E y 
and H z , is referred to as if-polarization. 

More specifically, as infers from the set of equations ffl3|) . in the case of 
-E-polarization the electric field is polarized along the ^-direction, while the 
magnetic field is confined to the plane XY. That is, E e>x = E ejV = H e z = 0, 
with the components of the magnetic field satisfying 

ie c 2 dE £tZ ie c 2 dE e>z 

u oy uj ox 

Substituting these expressions for H x and H y into the second line of ( TT3I) 
yields 

^ + ^ + ^ = »■ < 16 > 

E e = E e>z z, (17) 
H e = tf^x + H e , y y, (18) 

with E e>z satisfying Helmholtz's equation, according to ffT6]) . Analogously, in 
the case of i?-polarization the magnetic field is polarized along the z-direction 
and the electric one confined to the plane XY (i.e., H hx = H hy = E h>z = 0), 
with the components of the latter being 

„ i dHh z ^ i dHh z , . 

E h>x ~ , E h y - . (iyj 

ue Oy ue Q ox 



We thus have 



Substituting now these relations into the second line of f)T4j) yields 

8 2 H h . z d 2 H> 



M + + k 2 H h , z = 0, (20) 



dx 2 dy 2 



8 



which allows us to characterize if -polarization as 

Eh = Eh tX k + Eh >y y, 
Hh = H hz z, 



(21) 
(22) 



with Hh z satisfying the Helmholtz equation f )20|) . Therefore, any general 
(time-independent) solution will be expressible as 



E 
H 



E e + E h = E e + [V x H fc ] 

we 



H e + Hf, 



[V x E e ] + H 



h- 



(23) 
(24) 



Since and H z satisfy Helmholtz's equation, consider now that both are 
proportional to a scalar field, \I/(r) , which also satisfies this equation, i.e., 



E e = a\I/z, 

H fe = /3,/^e^z. 
V A*o 



(25) 
(26) 



Here, a and are real quantities and the phase shift between both compo- 
nents is given by 0; §5§ has been used to obtain the correct dimensionality 
in the r.h.s. of fl§6]). If (J25]) and ([26]) are substituted into Eqs. ([23} and ([1 
respectively, the latter become 



if3e irp d$> A iPe'tdV A 

E = — — — x — — y + a^z, 

k ay k ox 

ia A «a 9$ „ fc^e** TA 
H = -xH — y + — — ^z, 



with their time-dependent counterparts being 



(27) 
(28) 



E(r,t) 



H(r,t) 



x : — y + aWz 



ia 9* A 
xH — y 



-iuit 



ufj, dy Ufio dx 



^z 



-iuit 



(29) 



(30) 
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Equations ( ]2"9~]) and (13"U]) are general time-dependent solutions for a problem 
which can be described in terms of superpositions, as also happens with (1231 
and (j24p . Once this set of equations is set up, the whole problem reduces 
to finding \l/ and its propagation along x and ?/ (by the above hypothesis, 
the set of equations does not depend on z), which is a boundary condition 
problem. 



3. Incident EM plane wave, its polarization and photon trajecto- 
ries 

Before reaching the two slits, we assume that both the electric and mag- 
netic fields propagate along the y-direction. Moreover, we also assume that 
the incident scalar field is a monochromatic plane wave, 

*o(r) = e*». (31) 

Of course, this is just a simplification in order to develop some analytical 
expressions, where we have assumed that the incident wave is locally plane 
in a region that will cover the two slits. However, in a more realistic case, one 
could consider a more general incident wave function of a limited extension, 
like a beam or a wave packet. In the case of massive particles, the effects 
associated with incident wave functions of limited extension in two-slit ex- 



periments have been studied in ll|, |l2| (similarly, a study of the effects of 
the outgoing diffracted waves can be found in |10|). 
Introducing (ST} into ([29]) and flU} yields 

E (r,t) = + 

= [-fie l{ky+,p) k + ae lky z] e - *"*, (32) 
H (r,t) = [tfoA** + #o,e,*z]e- <u * 

— [ae iky it + Pe^+^z] e - *"*. (33) 




From these solutions, it follows that the polarization is going to play an 
important role in the interference patterns observed, and also in the topology 
displayed by the photon trajectories. 

Consider the electric field ( 132]) . whose real components are 

K Ax = -Pcoa(ky-ut + </>), (34) 
£J„ = a cos(% — tot). (35) 
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Since the magnetic field displays the same polarization properties as the 
electric field due to their relationship through the Maxwell equations (jl]) and 
(j3J) , we will only argue in terms of the electric field without loss of generality. 
Thus, expressing and (135]) as 



E r E r 

— lJL _| £:i cos — s in(fc?/ — ojt) sin0, (36) 

p a 

E r 

— — sin0 = cos(ky — ui) sin0, (37) 
and then squaring and rearranging terms, we reach 

' ~ \2 / — \2 / — \/~ 

Tpr \ ( rpr \ ( rpr \ ( Tpr 

-^OA^ \ I -^O.e.z 1 / -^0,^,2) 1 / "^0,6,2 \ , _ . 2 



According to this relation, several cases are possible depending on the value 
of the phase-shift 0: 

(a) When = or n, 

x 2 

| "^O.e^ \ _ „ _ ^0,h,x _ "^0,6,2 / Qn x 

— - — ± — u =^> — - — — =F • {.ov) 

pal pa 

This case describes linear polarization, for any arbitrary a and /3. 

(b) For any other value of 0, there is elliptic polarization. In this case, 
(138j) is the equation of an ellipse inscribed in a rectangle parallel to 
the XZ plane, with sides 2\a\ and 2\(3\ (see Fig. [1]). The electric field 
(and also the magnetic one, which is perpendicular) can move clock- 
wise or anticlockwise, as seen by an observer toward whom the EM 
wave is moving. These rotations define, respectively, right-handed or 
left-handed polarization states. The polarization handedness can be 
determined by computing the time derivative of the angle formed by 
the two components of the electric field, 9e = (tan) -1 (£?o, e ,z/-^oh x)> 

dO E 1 




dt 1 + (E^JE^ 

(3 sin cu [1 + tan 2 (% — ut)} 



a 1 + {%,JEI A 



) 2 

x) 



(40) 
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< > 

21/61 

Figure 1: Diagram showing the general elliptic motion described by the electric field due 
to its polarization (the same applies to the magnetic field, which is perpendicular). If 
4> = ±7r/2 and a — /3, the ellipse becomes a circle (circular polarization), and if <j) = or 
7r, it reduces to a segment (linear polarization) regardless of the a and /3 values. 



As can be noticed, the information about the handedness is contained 
in the prefactor of this expression, since the second factor (the time- 
dependent one) is always positive. Thus, if a and (3 are chosen positive, 
d0E/dt will be positive for — n < < and negative for < < ir. 
In the first case, the field is right-handed polarized (8e increases with 
time) and, in the second one, it is left-handed (9e decreases with time). 

(c) In the particular case <fi = ±vr/2, we have 

\ 2 

= 1. (41) 



rpr \ I rpr 



(3 I \ a 



If a = (3, the ellipse described by (jJTj) reduces to the equation of circle. 
We then have circular polarization. 

The particular values of 0, a and (3, which define the polarization state 
of the incident plane wave, are very important regarding the observation of 



interference patterns behind the slits according to the Arago-Fresnel laws (46 
But they are also going to be very important with respect to the topology of 
the corresponding photon trajectories, as shown below. 
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In Sec. IU we will tackle the calculation of the EM field and photon trajec- 
tories behind gratings with N and two slits. In the second case, we will not 
only consider free passage through the slits, but also when they are covered 
by linear polarizers parallel to the XZ plane, whose polarization axes are 
oriented along the z-axis in one of the slits and along the x-axis in the other. 
Moreover, and without loss of generality, we are going to assume that the po- 
larizer oriented along the z-axis produces E'-polarization, while the polarizer 
along the x-axis gives rise to if-polarization. As can be shown, the action of 
these polarizers on an incident plane wave which propagates along the y-axis, 
such as (l3~Tj) . justifies this specific labelling. Note that, after passing through 
such polarizers, the EM field described by Eqs. (13"2|) and ( 13"3"|) gives rise to 
transmitted fields, E tr and H tr , with the following polarizations: 

a) If the polarization axis is oriented along the x-axis: 

E tr (r) = -(3e^y + ^ = E 0Ax ic, (42) 

H tr (r) = .[^ (3e i ^z = H 0Az z. (43) 
V Mo 

b) If the polarization axis is oriented along the z-axis: 

E tr (r) = ae ikv z = E 0t<he z, (44) 

H*.(r) = .f^ae iky x = H 0te>x £. (45) 
V Mo 

Regarding the photon trajectories, before the EM field reaches the grat- 
ing, they are obtained from 

S (r) = l - Re[E (r) x H*(r)] 

= (« 2 + /? 2 )y, (46) 

E/o(r) = i[e E(r).E*(r)+/ioH(r).H*(r)] 

= |(« 2 + /? 2 ), (47) 

where the time-independent parts of ( |32|) and ( |33|) have been considered. 
Substituting ( T46l) and ( T471) into ( TTTT) . we obtain 

* = (48) 
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which, after integration, yields 



x(s) = x , z(s) = zq, (49) 
y(s) = y Q + s. (50) 

That is, the EME flow evolves along the y-direction and, therefore, photons 
pursue straight lines parallel to the y-axis. Since we are in vacuum, we can 
assume that the distance s travelled by a photon during a time t is given by 
s = ct and, therefore, ( 150]) can also be expressed as 

y(t) = y + ct. (51) 



4. EM field behind a grating 

According to Sec. [2J the EM field behind a grating can be described 
by Eqs. ( 127)) and ( 128)) (or Eqs. (129)) and (130)) . respectively, if we consider 
time-dependence), with the scalar function ^(r) = ty(x, y) satisfying both 
Helmholtz's equation and the boundary conditions at the grating. If the 
incident EM field is monochromatic, we can also assume that the incident 
scalar function is given by ()3ip . 

Traditionally, the exact solution for the slit array that we are going to 
consider here arises from Helmholtz's equation (see Eqs. (|6]) and (JTJ)) and is 



expressed Fresnel-Kirchhoff integral [49| 



= e- l/4 e ,ts I" + )e ik{x - x ' )2/2y dx' } (52) 

V 27T v J-™ 

where ^(x',0 + ) denotes the wave function right behind the two slits. Then, 
after assuming some appropriate approximations, expressions valid to de- 
scribe Fresnel and Fraunhofer diffraction can be obtained from the corre- 
sponding general solutions. Alternatively, Arsenovic et al. have shown 
that the solution behind a grating can also be expressed as a superposition 
of transverse modes of the fields multiplied by an exponential function of the 
longitudinal coordinate, i.e., 

= — L= e iky I c(k x ) e ikxX ~ ik * vl2k dk x , (53) 

V 27T J 



where 



c(k x ) = -= I V(x } + )e- tkxX dx. (54) 
v 2n 
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As also shown by Arsenovic et al. 5l|, [52), the solution (j53p is equivalent 
to the Fresnel-Kirchhoff integral (j52p whenever the wave number associated 
with the transverse mode k x of a general slit satisfies the condition k 3> k x , 
which occurs in most cases of physical interest. 

For a grating which is totally transparent inside the slits and completely 
absorbing outside them, the boundary conditions are: ^(x,0 + ) = for any 
x belonging to the slit support and ^(x,0 + ) = ^(x,0~) for any x within 
the apertures, with ^(x, 0~) being the wave function incident on the grating, 
here given by (13T|) . Thus, in the case of an incident beam illuminating iV 
openings of width S and mutual distance d, a simple analytical calculation 
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CN{k x ) 



2ttN 



sin(k x 5/2) 
kJ/2 



sm(Nk x d/2) 
sin(k x d/2) 



(55) 



In the particular case N = 2, i.e., the well-known double-slit experiment, it 
is useful to express Eqs. (1521) and (1511) as 



V{x,y) 



— e - i7T/4 e tky [ y(x',Q-)e ik{x - x ' )2/2y dx' 



V 2n v Ja 2 

fa(x,y)+ifa(x,y), 



(56) 



and 



C2{k x ) 



27C J a x 



y{x',0-)e~ ik * x dx 



V(x',0-)e- ik * x dx 



A 2 



~j= (c\,d/2 + Cl,-d/ 2 ) 



(57) 



respectively, where ?p\ refers to the scalar field coming from slit 1, centered 
at x = d/2, and ip2 is the scalar field coming from slit 2, at x — —d/2. After 
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carrying out each integral in ( 157)) . we obtain 

&m(k x 6/2) 



ci,d/ 2 (k x ) 



C\-d/2{k x ) 




2 

TlS 

2 

ix5 




sm(k x 8/2) 



A t 



-ik x d/2 



ik x d/2 



which, when they are added, yield 



c 2 (K 



sm(kJ/2) 
k x 8/2 



cos(k x d/2). 



(58) 



(59) 



In the space behind the grating, the EM field is given by Eqs. (f27|) and 
where \l/ is described by (15B|) . That is, the resulting EM field behind 
the grating consists of a superposition of two fields propagating from each 
slit, which reads as 



E 



k dy 
i/3e irt> dip 2 
k dy 
Ei + E 2 , 



k dx 
iPe** dip 2 
k dx 



y + a^iz 
y + aip 2 % 



(60) 
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ia 


dip! 




ia 


dipt 


y + 


kpe* 




dy 


x + 




dx 




ia 


dip 2 


x + 


ia 


dip 2 


y + 






dy 


<*>fJLo 


dx 





ipiz 
ip 2 z 



= Hi + H 2 . 
5. Photon trajectories behind the two slits 



(61) 



In order to obtain the photon trajectories behind the grating, first we 
substitute ( 127)) and (128]) into (jSj), which yields the components of the Poynting 
vector along the different directions, 



s z 



i(a 2 + (3 2 ) (_W 



dx dx 



— — 
dy dy 



4u;/i 

iap sin fdV d^* d^ d$>* 



(62) 
(63) 
(64) 
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Proceeding similarly with leads us to the EME density, 

4:U 2 fi \ dx dx dy dy J ' 

which describes the interference pattern at the observation screen. 

Before computing the photon trajectories, it is interesting to note the 
following feature about the interference pattern. Consider (165]) expressed in 
terms of the two scalar fields, ip\ and ifi 2 , i.e., 

(J = (a 2 + (3 2 ) fd^d^l | d^dr, | fe2 » 

4cu 2 /io V ® x dx dy dy 1 

(a 2 + /3 2 ) fdip 2 dip; dip 2 dip; 



Auj 2 /jo \ dx dx dy dy 
(a 2 + /3 2 ) fd^ dipt* d^ dip. 



4:U 2 fi \ dx dx dy dy 
{a 2 + (3 2 ) fdip 2 dip*i , dip 2 dipt 



+ a + V'Mi ■ ( 66 ) 



4oj 2 /i \ <9x cte ch/ dy 

In short-hand notation, (j6"Sj) can also be expressed as 

U = U 1 + U 2 + U 12 , (67) 

where JJ\ and U 2 are the EME densities associated with ip\ and ip 2 (the 
first and second terms in (166]) ). respectively. On the other hand, XJ\ 2 (the 
last two terms in (1661) ) is the EME density arising form the interference of 
these waves. Since the polarization part (prefactor in terms of a and /3) and 
the space part (depending on \1/) appear factorized in (165]) . the interference 
pattern observed will not depend on the polarization state of the incident 
field, in agreement with the Arago-Fresnel laws. 

Substituting now (j6"21) - (j6"5l) into ffTT]) renders the corresponding trajectory 
equations along each direction, 



dx 



Ts =ik \ d¥M~* X d^dv* X j2tt J> ((,(S) 

dx dx dy dy 

^ = ik I °y °i \ (69) 

ds ) Md^ + d*d^ + k2 ^ ? 

dx dx dy dy 
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d* d^* d* 



dz 2ial3sm.(j) I dx dy dy dx 



ds (a 2 + /3 2 ) 1 &B d** + d^d^ + fc2 ^. 



(70) 



dx dx dy dy 

As can be noticed from these expressions, all the information about the 
polarization state of the diffracted EM wave is contained in the prefactor 
of (170]) . Thus, regardless of the polarization of the initial EM wave, since 
the diffracted waves arising from each slit have the same polarization state, 
one will always observe interference fringes, which is in agreement with the 
Arago-Fresnel laws 



In the case of linear polarization, (1701) vanishes 45| and we can solve the 
photon-trajectory equations simply by parametrizing, for example, y as a 
function of x, i.e., 



d^* T *d^\ 
dy = V dy dy J 

dx ~ f^d^_ _ 

dx dx J 



(71) 



while the solution of ( I64p is simply z = z . On the contrary, in the case of 
elliptic polarization, the ^-component does play an important role, as can be 
noticed when the photon-trajectory equations are computed, 



dy V dy dy J 

dx f^^_y*^\ 
dx dx J 



d* d** d* dW 



dz 2a/3sin</> I dx dy dy dx 



dx (a 2 + /3 2 )k) ^d** 

dx dx 



(72) 



(73) 



In Fig. |2j the photon trajectories associated with the diffraction of an incident 
EM field circularly polarized (a = (3, = tt/2) are plotted. The parameters 
considered in the simulation are: A = 500 nm, d = 20A = 10 /zm and 
5 — d/2 — 5 /im. The projections of these flow lines on the XY plane, 
shown in Fig. 12(b), are identical to those for incident linearly polarized light 



45| (note that ( 1721) is exactly the same as ( 1711) ). As shown elsewhere 54 



within the context of Bohmian mechanics, for equations like (172]) . which 
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describes the photon-trajectory projections on the XY, trajectories exiting 
through one slit never cross the trajectories coming up from the other one. 
Moreover, if both slits are identical (here this means they have the same 
width and transmission function), the EME fluxes coming out from each 
slit are symmetric with respect to the axis y = 0. Now, although neither 
the electric nor the magnetic field depend on the z-coordinate, the photon 
trajectories display some remarkable features along this direction, as seen 
in Figs. [2]^c) and[2^d). This is an effect of having circular (or elliptical, in 
general) polarization, which vanishes in the case of linear polarization, when 
= or 7T. 

In order to understand the somewhat unexpected motion along the z- 
direction, let us go back to (16"4"I) . Rearranging terms and using (I6"2j) and (16"3"j) . 
this equation can be rewritten as 



a(3 sin< 



dS v dS, 



(a 2 + (5 2 )k \ dx dy 



a (3 sin < 



(a 2 + /3 2 )fc 



(74) 



where 



V 



x 
d_ 

dx 



y 
d_ 

dy 
S a 







(75) 



According to ( 1741) . the presence of a polarization state gives rise to a flow 
along the z-direction in terms of the vorticity manifested by the fields S x and 
S y , which may lead the photon trajectories to display loops out of the XY 
plane. Nodal structures and other sing ularities and topological structures 
may then appear, as shown by Nye 



55 



tions 



56 



within the context of wave disloca- 
and by other authors within the context of the Riemann-Silberstein 

see Appendix IA~j) . 



38. 57. 58 



complex formulation of Maxwell's equations 
Experimentally, what one would observe on the XZ plane is simply the 
typical fringe-like interference pattern constituted by dark and light paral- 
lel strips, which results from the accumulation of photons arriving at this 
plane. Note that ( 165]) describes the interference pattern, which results from 
transporting the EME density through space (from the slits to some detec- 
tion screen) in accordance to the guidance or continuity equation ( ITOl) . This 
means that, if we make a histogram with the arrivals of a statistical distribu- 
tion of photon trajectories along the x-direction (from now on, we will label 
the normalized amount of such counts as £,-, with i — x, z), the well-known 
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Figure 2: Photon trajectories (15 for each slit) behind a two-slit grating associated with 
an incident EM plane wave is circularly polarized (a — /3, <p — 7r/2): (a) 3D view, (b) XY 
projection, (c) XZ projection and (d) YZ projection. The parameters considered in the 
simulation are: A = 500 nm, d = 20A = 10 fim and S = d/2 = 5 /im. 
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interference pattern emerges, as can be seen in Fig. [3](a). However, from (175 p . 
all the arrivals at a certain height zt will not arise from positions at the slits 
at the same height zq = Zf, but there is a flux upwards and downwards which 
breaks the longitudinal (along the ^-direction) symmetry of the experiment 
when it is studied from the viewpoint of photon trajectories. This gives rise to 
a certain distribution of arrivals along the z-direction, as shown in Fig. E^b). 
Since the photon trajectories distribute evenly around z (here, along posi- 
tive and negative z, since we have chosen zq — 0), as can be appreciated in 
Fig. E^d), their distribution is also going to be symmetric with respect to z 
in the histogram of Fig. [3]^b). Nonetheless, we would like to point out that 
this effect along the z-axis arises because we are considering a particular ini- 
tial value of the z-coordinate. If this coordinate is sampled properly up and 
down, this effect cannot be appreciated, but only a homogenous distribution, 
unlike what happens along the x-coordinate. 

6. The EM field and photon trajectories behind two slits, each 
followed by a linear polarizer 

According to the Arago-Fresnel laws [46|, if two diffracted beams with 
different polarization states interfere, the visibility of the interference pattern 
will decrease. Indeed, if the polarization states are orthogonal the pattern will 
disappear totally. In order to describe this effect with photon trajectories, 
we consider the EM field behind a grating with two slits, each followed by 
a linear polarizer, such that behind the slit 1 there is a polarizer with its 
polarization axis oriented along the z-axis and behind the slit 2 the polarizer 
is oriented along the x-axis. As in Sec. HI we will also start expressing the 
total EM field in terms of the superposition of two fields, each propagating 
from one of the slits (see (160]) and fl6~Tj) ). However, due to the presence of 
the polarizers and their filtering effect on the incident EM field, instead of 
Eqs. f l2l)|) and f l2"E|) . now we will express the electric and magnetic fields in 
terms of ipi and ip2, described by (!56|) . i.e., 



E e = a^i(r)z, 



(76) 




(77) 



21 



(a) 




Figure 3: Histograms built up by counting the end points of individual photon trajectories 
associated with the incident EM field circularly polarized (a = /3, <fi = tt/2) of Fig. [51 The 
detection screen is at L = 1 mm from a two-slit grating and the histograms represent 
counts: (a) along the x-direction and (b) along the z-direction. In part (a), the red solid 
line indicates the theoretical curve predicted by standard Electromagnetism, according to 
dHSJ). The parameters considered in the simulation are: A = 500 nm, d — 20 A = 10 /xm 
and 5 = d/2 = 5 /im. The total number of trajectories considered is 5,000, with initial 
conditions homogenously distributed along each slit. 
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By substituting <^ and (JTZD into (J23D and (JMD, we will obtain the EM field 
behind the two slits covered by the polarizers, 



E = — — x — y + afaz, 78 

k dy k dx 

ia dipt ia d^ kpe^ A 
H = — xH — yH ^2Z. (79) 

From these fields, we obtain the expressions for the EME density and the 
Poynting vector, 

4:U 2 fi \ dx dx dy dy 1 



, /3 2 <¥| + 9^ + fc2 , 



tfi + C/ 2 (80) 



and 



ia 2 / <9^* „ dipi \ i/3 2 f dip2 * dip 



4a;/io V 1 dx J 4cl>/z \ dx 2 dx J ' 

ia 2 ( Ml di;\ if? ( M>1 ,Mh\ .... 

Aku^io V <9x <9x 9y / 

AkuJiiQ \ dy dx dx dy J ' 

respectively. As can be noticed in ( |80|) . £/ describes the bare addition of EME 
densities associated with the fields diffracted by each slit, with no interference 
term. This means that no interference pattern is going to be observed at the 
detection screen, in accordance to the Arago-Fresnel law for the interference 
of two beams with orthogonal polarization states (perpendicularly polarized 
in the case of linear polarization |47f or with opposite handedness in the case 
of elliptic or circular polarization [48]). Regarding the Poynting vector, we 
note that the EME flux along the x and ^/-direction is also given by the simple 
addition of fluxes coming from each slit. Thus, unlike the diffraction problem 
dealt with in the previous section, now it is not possible to factorize the flux 
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(neither the EME density) in terms of its spatial and polarization parts. This 
has a consequence on the EME flux along the ^-direction, which does not 
satisfy the rotational character described by (I75p . On the other hand, since 
the third and forth terms in fl83|) do not vanish, there will be an EME flux 
along the z-direction even in the case of linear polarization (remember from 
Sec. E]that, for linear polarization, S z = 0). This can also be seen from the 
photon-trajectory equations, 

iaf3e l<t> ( dip 2 dtp{ dip 2 dtp\ 



dz k \ dx dy dy dx 



dx 2 / dipl ^dipA . r2 (, di)* 2 dip 2 



dy dy J \ dy dy 

ia/3e~ l,l> f dipi dipl dipt dip 2 



k \ dx dy dy dx 



1 dy 1 dy J \ 2 dy 2 dy 



(84) 



2 / , ,*drk\ 2 / d^ 2 d^j 2 

_ \ ft; dy J y dy dy 

a ^"5 V'l-S- +Z 3 ^"5 ^2 



<9x <9x / V dx dx 



(85) 



Note in (}84"|) that dz/dx, effectively, does not vanish, not only for <fi = or tt 
(the conditions for linear polarization), but neither for any other value. 

In Fig. HI we observe an ensemble of photon trajectories for an incident 
EM field linearly polarized (a = ft, (j) = 0). As before, the parameters 
considered in the simulation are: A = 500 nm, d = 20A = 10 /im and 5 = 
d/2 — 5 /im. As can be noticed, when polarizers with orthogonal polarization 
directions act on the diffracted wave, the topology of the flow lines changes 
dramatically when compared with that observed in Fig. [2J When one looks 
at the XY projection (see Fig. 11(b)), the wiggling behavior that gives rise 
to the different interference fringes of the pattern in the Fraunhofer region 
are lacking. This reflects the fact that the EME density, described by (180|) . 
is just the sum of the EME densities associated with the components of the 
wave that arise from each slit, which can be appreciated in the histogram 
presented in Fig. [5(a). This histogram reproduces the distribution pattern 
of the photon trajectories in the Fraunhofer region, which simply consists of 
the direct sum of the two single-slit diffraction patterns associated with the 
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Figure 4: Photon trajectories (15 for each slit) behind a two-slit grating, where each 
slit is followed by polarizers with orthogonal axes and the incident EM field is linearly 
polarized (a = 0, <fi = 0): (a) 3D view, (b) XY projection, (c) XZ projection and (d) YZ 
projection. Because the incident EM field is linearly polarized, the ensembles leaving each 
slit behave exactly the same and, therefore, the trajectories exiting through slit 1 look the 
same as those exiting through slit 2. The parameters considered in the simulation are: 
A = 500 nm, d = 20A = 10 ^m and S = d/2^ 5 fim. 



EM field arising from each slit. In particular, since the grating is assumed as 
fully transparent along the space covered by the slits and opaque elsewhere, 
the single-slit diffraction pattern is a square sine function, this leading to the 
appearance of the maxima and minima observed in Fig. [5(a). The reason 
why this pattern looks like the one arising from a single slit rather than the 
sum of two of them is very simple. The observation distance L (1 mm) is 
much larger than the distance between the two slits d (10 /im) and, therefore, 
the single-slit diffraction pattern associated with each slit basically overlap. 
For example, the first minima for slit 1 occur at x]^[ Q ~ —105 /im and ~ 
95 /zm, while for slit 2 they are at x^ in ~ —95 /im and x^ in ^105 pm; 
therefore, the minima of the total pattern will appear at x^ in ~ ±100 /im, 
as can be seen in Fig. 0(a). It is worth mentioning that the loss of the 
interference fringes as well as the lack of wiggling features in the trajectories 
is similar to the behavior found in trajectories for massive particles in the 



case of decoherence 59, 60 



If we look at the topology of the photon trajectories along the z-direction 
(see Figs. H(c) and (d)), we note that there is a symmetry with respect to 
x = (see Fig. 11(c)). This is the result of the breaking of the rotationality 
property mentioned in the previous section. The same "symmetry" breaking 
can also be observed by looking at the distribution of the trajectories along 
the z-direction, as shown in Fig. [5(b). If we changed the initial polarization 
state, from = to = tt, we would obtain the same pattern, but inverted 
with respect to z = 0. Again, as we discussed in Sec. this is an effect 
arising from considering a particular initial value of the z-coordinate for the 
photon trajectories, which will disappear when a sampling along the z-axis 
is considered. 

In the case of an incident EM field circularly polarized, displayed in 
Fig. El we find that the photon-trajectory projections on the XY plane (see 
Fig. [6(b)) look exactly the same as those for the incident EM field linearly po- 
larized (see Fig. E(b)), in accordance to Eq. (1851) . However, as seen in Sec. [5j 
due to the rotationality associated with circular polarization, a breaking of 
the specular symmetry with respect to x = will take place, and the EME 
flux leaving each slit is going to be different, as shown in Figs. [5(c) andE(d). 
This will give rise to a histogram like the one shown in Fig. [5(a) when pho- 
ton are counted along the x-direction and another similar to that of Fig. [3(b) 
when they are counted along the z-direction. 
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Figure 5: Histograms built up by counting the end points of individual photon trajectories 
associated with the incident EM field linearly polarized (a = /3, <fi = 0) and where, as in 
Fig. |4l each slit is followed by polarizers with orthogonal axes. The detection screen is 
at L = 1 mm from a two-slit grating and the histograms represent counts: (a) along 
the x-direction and (b) along the ^-direction. In part (a), the red solid line indicates 
the theoretical curve predicted by standard Electromagnetism, according to ([80jl . The 
parameters considered in the simulation are: A = 500 nm, d — 20A = 10 /mi and 5 = 
d/2 = 5 fim. The total number of trajectories considered is 5,000, with initial conditions 
homogenously distributed along each slit. 
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Figure 6: Photon trajectories (15 for each slit) behind a two-slit grating, where each 
slit is followed by polarizers with orthogonal axes and the incident EM field is circularly 
polarized (a — (3, <p — k/2): (a) 3D view, (b) XY projection, (c) XZ projection and (d) 
Y Z projection. The parameters considered in the simulation are: A = 500 nm, d = 20A = 
10 ixm and S = d/2 = 5 /im. 
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7. Final discussion and conclusions 



Interference experiments with non-relativistic massive particles and pho- 
tons can be well understood on the theoretical grounds provided by the (non- 
relativistic) Schrodinger and Maxwell equations, respectively. Thus, if one 
wishes to determine photon trajectories on the same footing as Bohmian 
trajectories, the most appropriate theoretical framework is the one based 
on Maxwell's equations, more specifically their hydrodynamical formulation, 



proposed by Bialynicki-Birula |27l . |28| . |29| . |30|. Within this formulation 



close analogy can be established between the trajectory equations based on 
the latter and Bohm's approach when the particle aspect of both light and 
matter is taken into consideration. This allows to compare on the same 
grounds the (Bohmian) trajectories for massive particles with the trajecto- 
ries derived for photons from classical electromagnetism. The photon trajec- 
tories (EME flow lines) are determined from the Poynting vector, with the 
components of the electric and magnetic vector fields expressed in terms of a 
function that explicitly takes into account the boundary conditions imposed 
by the grating. It is remarkable that, in the case of photons, any trajectory- 
based interpretation will be complementary to the standard Huyghens' one, 
based on the superposition of secondary wavelets. Furthermore, the topology 
displayed by the photon trajectories is strikingly similar to that displayed by 
massive particles. As happens in the case of massive particles |l4|, such a 
topology can also be inferred and explained from the corresponding trajec- 
tory equation. 

Here, we have considered the Bialynicki-Birula hydrodynamical frame- 
work to analyze the effects of polarization on interference pattern for waves 
which arise from an initially polarized monochromatic EM plane wave. Ac- 
cording to the Arago-Fresnel laws of interference with polarized light in 



Young-type experiments |46|, if, for example, the two interfering beams are 



linearly polarized in orthogonal directions 47| or both are elliptically po 



larized but with opposite handedness |48|, no interference will be observed 



The photon trajectories presented in this work show how the EME density 
distributes in configuration space, giving rise to the appearance or not of 
interference fringes. This thus constitutes a very interesting corpuscular-like 
description of the Arago-Fresnel laws for the interference of polarized light, 
consistent with the more standard wave description. Moreover, it is clear 
that the photon wave function can be interpreted as an energy localization 
amplitude or detection amplitude within the hydrodynamical approach to 
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Electromagnet ism. Now, from the agreement between the histograms and the 
corresponding EME density distribution (which follows from fllOp or fill]) ), 
this density can also be understood as a probability density function in the 
sense that arrivals are going to distribute statistically according to it, as can 
be seen in the recent experiment carried out by Dimitrova and Weis |4], for 
example. 
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A. Appendix A: Photon trajectories within the Riemann-Silberstein 
formulation of Electromagnetism 

Both Maxwell's and Schrodinger's equations describe the evolution of 
fields in configuration space, namely the EM field and probability fields asso- 
ciated with non-relativistic massive particles, respectively. From these fields, 
one can obtain trajectories which show how they evolve in space, namely pho- 
ton trajectories or Bohmian trajectories, respectively. However, Maxwell's 
equations look quite different formally from the Schrodinger equation as they 
are usually given. A formulation that allows one to put Maxwell's equations 
on similar formal grounds as Schrodinger's one (but without H due to the lack 
of mass of photons) is the complex form of Maxwell's equations 43j,|44|, which 
is based on the so-called Riemann-Silberstein complex EM vector 




F(r,t) 



V2 



e E(r,t) + i v ^H(r,t) 



(86) 



with the change of variable 



E(r,t) 
H(r,t) 



1 



iy/2jIo 



F + F 
F-F 



(87) 
(88) 
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where E and H are real fields. Introducing ( 18 7p and ( 18 8 p into Maxwell's 
equations in the absence of electrical charge densities, we obtain 

i — = cV x F (89) 
V-F = 0. (90) 

As can be noticed, (1591) is the analog for photons of the Dirac equation for 
massive particles, while (I9"0"j) describes the conservation of the EME density 
through space. This analogy within the Riemann-Silberstein formulation 
becomes more apparent by gathering (|89|) and (|90|) in a single equation. 
This is done by applying the operator —id/dt to both sides of (|89|) and then 
rearranging terms taking into account ( 1901) and the vectorial relation 

V x (V x A) = V(V ■ A) - V 2 A, (91) 

where A is a general vector field. This renders 

d 2 F 



dt 2 



c 2 V 2 F, (92) 



which has the well-known form of the Klein-Gordon equation. In the particu- 
lar case of diffraction problems, which can be reduced to boundary condition 
problems because of their time-independence, the space and time parts of 
(|92|) are separable. Thus, Fj can be decomposed as 

Fi(r, *) = F<(r) = F<(r) e^, (93) 

where the space part, Fj, satisfies Helmholtz's equation, 

V 2 Fi(r) + k 2 Fi(r) = 0, (94) 

and its time-dependent part the differential equation 

d 2 Mt) 



dt 2 



-oj'cPM, (95) 



with c — u/k. 

Within this formulation, the EME density is given by 



U = - (e E ■ E + hq H • h) = F ■ F* (96) 
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and its flux, described by the Poynting vector, is straightforwardly obtained 
after developing the time-derivative of U to yield 



S = E x H = icF x F* 



(97) 



where the relation 



V(A xB) = B- VxA-A-VxB 



(98) 



for any two vectors, A and B, has been taken into account. From these 
magnitudes, we can now obtain the stationary photon trajectories within 
the Riemann-Silberstein formulation as, 



(here, (A) denotes the time-average of the magnitude A), which transport 
the time-averaged EME density, U, as described by the (also time-averaged) 
Poynting vector, S. 

If the electric and magnetic fields are complex and the definition of the 
Riemann-Silberstein vector is kept as in (I86p . i.e., the real and imaginary 
parts are given by the electric and magnetic fields, respectively, then we 
need to include into the formulation two of these vectors in order to have a 
complete description of the problem, each one associated with the real or the 
imaginary parts of the electric and magnetic fields. That is, if 



with Ej and H,; (i = 1, 2) being real vector fields satisfying the corresponding 
Maxwell equations, we will have 



dr _ 1 S _ .(F x F*) 
ds~ cU ~ (F- F*) 



(99) 



E = Ei + «E 2 



H = H x + m 2 



(100) 




(101) 



The EME density (f96|) and the Poynting vector (1971) then read as 



U = ~(e E-E* + /i H-H*) = ]TBVF*, 



(102) 



i=l,2 




(103) 



i=l,2 
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respectively, and their time-averaged homologous, assuming the decomposi- 
tion ( I9"3"j) . as 

i=l,2 

S = ~ Re ^ic F i x F * J • ( 105 ) 

Apart from the interest of this formulation within the field of the Fun- 
damental Physics, it has also been considered in a more applied way, for 



example, in Solid State Physics and Condensed Matter [61|, 162 
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